preprint: NSF-KITP-07-01 



On the Sensitivity of a Hollow Sphere as a 
Multi-modal Resonant Gravitational Wave Detector 



F. Dubath^ J. Extermann^ and L. GottardP 

1 KITP, Kohn Hall, UCSB CA 93106, USA 

^ GAP - Biophotonics, Universite de Geneve, 20 rue de I'Ecole de Medecine, CH-1211 
Geneve 4, Switzerland 

■^SRON National Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, 
Netherlands 

E-mail: dubath@kitp.ucsb.edu 

Abstract. We present a numerical analysis to simulate the response of a spherical 
resonant gravitational wave detector and to compute its sensitivity. Under the assump- 
tion of optimal filtering, wc work out the sensitivity curve for a sphere first taking into 
account only a single transducer, and then using a coherent analysis of the whole set 
of transducers. 

We use our model for computing the sensitivity and therefore compare different 
designs of spherical detectors. In particular we present the case of 1 meter radius bulk 
and hollow spheres equipped with transducers in TIGA configuration, and wc explore 
the sensitivity of a hollow sphere as a multi- modal detector. 



PACS numbers: 04.80.Nn, 95.55.Ym 



Submitted to: Class. Quantum Grav. 



X On leave from DPT, Universite de Geneve, Switzerland 



On the Sensitivity of a Hollow Sphere as a Multi-modal Resonant Gravitational Wave Detector^ 
1. Introduction 

Direct detection of gravitational waves (GW) is still an unachieved goal. The sensitivity 
improvement of single detectors [1], to achieve the first detection, goes in parallel 
with the set-up of detector networks to perform real GW astronomy [2]. To that 
purpose, future GW detectors should be able to track GW arrival directions. Due 
to their symmetry, spherical resonant detectors are natural GW 'telescope' candidates. 
Moreover, such detectors are based on the same technology as resonant bars and can 
benefit of the decades of experience gained by the GW community with this kind of 
experimental setup. 

A technical problem resides in the fact that big cryogenic resonant spheres with 
high quality factor are difficult to build and to cool down. One can partially bypass this 
problem using hollow spheres and we will eventually focus on this kind of resonator. 

However, the main limitation when dealing with resonant detectors, is the small 
bandwidth compared with laser interferometers [3]. 

Great progress have been made using the tuning of the electrical mode [1] which 
allows to enlarge the bandwidth from few Hz to the order of 100 Hz around the resonance 
frequency. A possible way to further enlarge the bandwidth consists of monitoring 
other resonance frequencies of the antenna. In the case of a cylindrical detector this 
strategy is not convenient since the coupling of the n-th longitudinal mode to GW falls 
off as (2n + However for a sphere the next mode quintuplet is tightly coupled 

to GW. Monitoring these modes is possible using many transducers, tuned at different 
frequencies. Here we explore this possibility and also consider a model of a double-mode 
transducer sensitive to two modes at different frequencies. We have to stress that both 
issues are experimental challenges but they can lead to great improvement for spherical 
detectors. 

It is possible to enlarge the detector bandwidth by monitoring also the second 
quadrupolar mode. In this way one could for example study the radiation emitted by 
a binary system, consisting of either neutron stars or black holes, in the in-spiral phase 
and determines the chirp mass by measuring the time delay between excitations of the 
first and second quadrupole modes of an hollow sphere [5]. Used as a multi-modal 
detector even a single sphere is able to set limit on the stochastic GW background 
[7]. Furthermore, with an appropriate resizing of the sphere, the second quadrupolar 
modes can be shifted to the frequency region where existing small spherical detectors 
are sensitive. This would open the possibility of coincidence search between several 
spherical detectors and the DUAL detector [21], building in this way the base for a 
powerful omnidirectional gravitational wave observatory. 

In order to be able to compare different GW detectors we need to work out the 
sensitivity with a "standard" approach and we decided to use the strain sensitivity. 
A first goal of this paper is to furnish a mathematical framework to describe all the 
internal noises of a given detector and its response to an excitation (as example to a 
GW burst). This is realized as a matrix model in the Fourier space. Using this model, we 
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set up numerical simulations and first compute the sensitivity of each single transducer 
attached to the sphere. We are fully aware that performing the data analysis in this 
way does not exploit all the capacity of the spherical geometry. However, this first 
computation is instructive and may be useful in order to calibrate a future experiment. 
We point out the fact that the noise coming from the other transducers will then be 
an important limitation to the sensitivity. A clever way to use a spherical detector is 
to perform the data analysis combining the signals from all the transducers. Such a 
strategy, which requires multidimensional matched filtering [IZl [HI [I3] , exploits all the 
geometrical properties of the detector and therefore gives best sensitivity. This is the 
sensitivity to be compared with other existing computations [H]. 

The second part of this paper will be dedicated to the application of our model 
to "realistic" detectors having in mind the comparison of different possible spherical 
resonant GW detectors. We will focus on hollow spheres and in particular we explore 
the different way to realize a multi-modal detector. The possibilities we test are all 
technical challenges (dealing with 12 transducers, transducer inside the hollow sphere, 
double mode transducers) but we eventually show that there is no gain in working with 
transducer inside the sphere and that 6 double mode transducers may compete with 12 
single transducers configurations only if they have a very peculiar design. 

2. First part: the model 

2.1. Description of spherical resonant GW detectors 

2.1.1. The modes of a sphere The vibrational motion of a rigid body can be split 
into eigen-modes. For a sphere there are two families of vibrational eigen-modes: 
toroidal and spheroidal [in]. Each mode can be described as a forced damped harmonic 
oscillator, which means that the j-th mode [§| amplitude Zj satisfies in Fourier space the 
equation (see for example [6]) 



where is the physical mass of the sphere, ijJs,j is the j-th mode eigen-frequency, Qsj 
its quality factor and Fj{uj) the Fourier component of the forcing. The forcing of the 
mode is due to a stochastic (Langevin) force and to external forces. The stochastic force 
corresponds to the thermal excitation of the sphere mode. Among the external forces 
we retain only the tidal forces induced by GW, and forces due to the coupling of the 
mode with the transducers. 

We use the frequencies ujs,j as an input for our model. The computation of the uOsj 
can be found in [TT] . 

§ We use j to describe collectively all the numbers needed to specify the mode. For spheroidal modes 
j={n, m} [5]. In the following we will work with a subset of the modes. In particular we will use only 
J spheroidal modes that we will label by j = 1, .., J 




(1) 
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Figure 1. The transducer model. 

2.1.2. Transducer To give a clear presentation of the noises and the sensitivity we need 
a full model for the transducer. We use capacitive transducers (Figured]), based on the 
one employed on bar detector [1] and the spherical detector MiniGRAIL [12], composed 
by one resonator coupled to a dc-SQUID with input transformer, where the SQUID is 
described as a linear current amplifier [1^. However our formalism can be adapted for 
other kinds of transducer. 

Each transducer (a mechanical resonator with its read-out) will be modelled by a 
set of p = 3 coupled differential equations, driven by the specific intrinsic noises and its 
coupling to the sphere. In Fourier space each transducer and each read-out are described 
by the equation set 

mt{uj'i-uj^ + i'^) -lE/u //\ 

E r + i{ujL^-j^) -iMuj / I = K (2) 

-iMuj iuj{Ls + Li)) \ij \Vn) 

where m^, ujt, Qt are the resonators mass, eigen-frequency and quality factor, E is the 
electric field in the capacitor Ct formed with the resonator, r the transformer resistance, 
Li, Lp, Lg the inductances of the SQUID input coil, primary and secondary of the 
transformer, and M the mutual inductance of the transformer; see Figure [H x is the 
transducers position [T5|, Ip, li are the currents into the transformer and the SQUID, 
and li is the measured quantity. Finally, Vr and Vn are the voltages corresponding to the 
intrinsic noise and / is the sum of the thermal force and the forces due to the coupling 
to the spheres modes, which are described in Section 12.1.41 

In what follow, we will use a set of transducers and therefore add an index 
k = 1, .., N to all the quantities describing the transducer. 

2.1.3. Double-mode transducer Having in mind the possibility of building a multi- 
modal spherical resonant GW antenna we note that a simple way to achieve this goal 
is the use of multi-modal transducers. A multi-modal transducer has to be designed in 
order to be resonantly coupled to two distinct frequencies. As a model, one can think 
to a multi-modal resonator with different vibration frequencies, (each frequency tuned 
on one of the sphere resonance frequencies) and corresponding electronic read-out. The 
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resonator can be modelled as two oscillators located at the same point of the sphere. 

In Figure [21 we show the schematic diagram of a capacitive double mode transducer 
composed of a resonator mechanically coupled to the two first spheroidal modes of the 
sphere and a single SQUID read-out. The two resonating mass schematically drawn in 
the figure should be considered fixed at the same point on the sphere surface. Such a 
transducer could be fabricated using, for example, a mushroom type resonator placed 
inside a ring-shaped membrane transducer or using the geometry suggested in [25] . 
By means of two electrically insulated electrodes and two super-conducting matching 
transformers, one is able to couple the signal from both the resonators to a single 
SQUID amplifier, simplifying greatly the detector read-out electronics and the cryostat 
wiring. To achieve an optimal impedance matching each electrical mode defined by the 
two LC circuits should be coupled to the two corresponding spheroidal modes. In this 
scheme a mixing may occur between the two modes. However we will show here that 
this will not significantly reduce the detector sensitivity. 

The double-mode transducer is described by a set of p = 5 equations: two for the 
mechanical modes, two for the transformer currents /i, I2, and one for the SQUID input 
current 
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2.1.4. Coupling transducers to the modes of the sphere Since none of the state-of-the- 
art transducers are able to couple efficiently to tangential motion, we will restrict ourself 
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Figure 2. The double transducer model. 

to transducers only sensitive to radial motion of the sphere, and consequently we will 
neglect the coupling to toroidal modes [10] . 

We call {9k, (pk) the location of the k-th. transducer. At this position on the sphere 
surface, the j = {n, /, m}-th spheroidal mode induces a radial displacement described 

ani{Rs)Yim{Ok, (pk) = ajBjk (9) 

where aj = ani{Rs) is the radial eigen-f unction evaluated at the sphere surface, and 
is a spherical harmonic. This equation defines the pattern matrix Bj^. Using this 
matrix we obtain the equation of motion ( [6], [I3]). The coupling between the sphere 
modes and the transducers is then given (in the Fourier space) by 

= Fr'^'''^ + E c^^B.k ( {u;l + m,,5. - /—) (10) 

fk = fr^' + m,kuj' E B.kafz, . (11) 
j 

We obtain a set of {J + pN) coupled equations, with J the number of modes taken 
into account that we discuss bellow, the number of transducers and p the number of 
equations need to describe a single transducer. 

It is important to note that generally a transducer is coupled with a subset of 
the five sphere modes and therefore the presence of transducers provides an indirect 
coupling between the different modes. Reciprocally, different transducers are coupled 
through the sphere modes. Therefore, the total noise spectrum of a given transducer 
has a contribution from the intrinsic noises of the other transducers. 

We are principally interested in the quadrupolar {£ = 2) spheroidal mode family. 
The lowest quadrupolar multiplet {n = 1,£ = 2) has the advantage to contain the 5 
spheroidal modes with the lowest resonance frequency [11] . When we are only interested 
in those modes we can reduce the sphere vibration to these 5 modes and neglect the 
effect of the modes at higher frequency. When we describe a multi-modal detector, 
we are also interested in the second quadrupolar multiplet (n = 2,£ = 2) at about 2 
times the fundamental mode resonance. We note that other modes, not coupled to GW, 
{n = 1,£ = 0, 1,3,4) have their resonance frequency between the ones of the first and 
second quadrupolar multiplet [H]. These modes induce additional thermal noise sources 
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and have to be included when calculating the sensitivity of a multi-modal resonator, like 
previously shown in [22] . 



2.1.5. Bulk sphere model The remaining of the computation is easier to follow if applied 
to a concrete example. For this purpose, we present a simple model with only the 
spheroidal quadrupole modes of the sphere (J = 5, 1 = 2), and with = 6 capacitive 
transducers placed in the TIGA configuration (JJ] [6] . Note that we can write the system 
of ( J + pN) equations in the same form independently of the number of modes and 
transducers. We also use this example in order to set the notations. 

Collecting the equations f|T][2][9][T0|) and ffTTl) we obtain a description of the entire 
detector as 
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where z = 






: , and so on for the other variables. 
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S is a 5 X 5 diagonal sub-matrix given by the left hand side (LHS) of([T]) ,Tisa3-6x3-6 
sub-matrix with structure given by the LHS of ([2]) (each number in ([2D is now a 6 x 6 
diagonal matrix). The matrices Ci and C2 are 5x6 (resp. 6x5) sub-matrix read out 
of fllO|lll) which describe the mechanical coupling between the spheroidal modes and 
the transducer modes and are given by 
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2.2. Detector noise description 
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Starting from equation (1121) . knowing the forces and the voltage acting on the detector 
allows us to compute the sphere modes z, the displacement q of the transducer, and the 
currents Ip and Ij. To do that we invert the Z matrix and rewrite f[T^ as 

/ Z \ / pnoisG+GW \ 

^noisc 



q 




(15) 



II The precise location of the transducer changes the B matrix causing to loose the TIGA configuration. 
However this do not change the following analysis 

% As we have only modes with I = 2, all the aj are equal. Therefore, we drop the indices. 
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Mode amplitudes, transducer displacement and transformer currents are not directly 
measured and we are therefore interested only in the N last lines of G, which give the 
proportionality coefficients between forces (and voltages) and SQUID input currents Ij. 
We now describe the different noise contributions. 



2.2.1. Noise description We restrict ourself to the case where the disturbances 
pnoisc^ poise are only due to thermal excitations. In this case we can only 

access statistical property of these forces that is: 

= iy{t)) = , (16) 

{F{t)F{t')) = Ao6{t - t') {V{t)V{t')) = Al6{t - t') . (17) 

Furthermore one can compute [I6] that the coefficient Aq takes the form 

= 2kBTm/Q A^ = 2kBTr . (18) 

where is the Boltzmann constant, T is the thermodynamic temperature of the 
detector, and r the transformer resistance. 

Using the definition of the single-sided spectral density of the force which is obtained 
through the autocorrelation 

{Fimt')) = I r |^5H^)e-^-(*-*') , (19) 
Z Jo Ztt 

we obtain 

ksT 

Spj = ArUsUJsj— — (20) 

^ s,j 

S/,fc = ^rnt^k^t,k^^ (21) 
SK,fc = AkeTru. (22) 

The SQUID has also intrinsic noises which can be split into voltage and current 
noises. The determination of these noises requires the knowledge of the complete SQUID 
design. Therefore, we have to go beyond the description of Section 12.1.21 In particular 
we need the shunt resistance Rsh,k used to remove hysteresis, the washer inductance 
LsQ,k, and Msq^k the mutual inductance between the SQUID input and the washer, see 
Figure [3l 



9 




Figure 3. SQUID detail. 
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The SQUID voltage and current noises are, according to Clarke's model [H] 

Sv„,. = ll^f^-^Ml^,, (23) 

J^sh,k 

' LsQ,k \ kpTsQ^k 

^SQ,k J Rsh,k 

where Tsq^k is the SQUID thermodynamic temperature. 



= 16 ■"^--^'^ (24) 

\MsQ,k / Rsh,k 



{u) + iwiuj) = Gi{u)F{u) + iw{uj) (25) 



2.2.2. Noise matrix The noise transducer outputs are proportional to the SQUID input 
current and thus given by 

/ jpnoisc \ 
£noise 

V v„ ; 

where Iw{^^) is the SQUID noise current and Mj is a iV x (37V + J) matrix given by 
M, = (OAr,j On,n On,n ^n) (26) 

used in order to conserve only the SQUID input current and we have defined Gj{uj) = 
MjG(ci;) to simplify the notation. 

Knowing the current we can compute the noises matrix S , which is given by 

S = li(^)IiH^) = G/FFtGt + GiFll + iwF^G\ + iwi\v (27) 

The matrix FF"!" is the correlation matrix of the different forces and voltages (Note 
that not all entries have the same units). Assuming that the different forces and voltage 
are only due to noises and are not-correlated this reduce to a diagonal matrix containing 
only autocorrelation 

FFt = Diag(SF,Sf,Sv„SvJ (28) 

The same argument leads to the cancellation of the Fi|y and IwF^ terms into ( 1271) and 
to set 

iwil = Diag(SH.) (29) 

Therefore under the assumption that the different noises are not correlated the 
noise matrix (1271) reduces to 

gnoise ^ G.Diag (Sp, Sf , Sv., SvJ Gl + Diag(SH.) (30) 
2.3. Effect ofaGW 

The presence of a GW will manifest itself as a force acting on the sphere modes. We 
skip the computation (see [6]) and just note that in Fourier space the force acting on 
the j-th mode is 

w = A^^rrisXjRshj (31) 
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where XjRs, the effective length of the mode, depends only on the multiplet (n) to which 
the mode belongs, and hj is the projection of the (spatial part of the) GW tensor El on 
the j-th mode. Choosing the decomposition of the tensor on real matrix [15] we specify 
the form of h for the quadrupolar modes {i = 2) b}0. 



V3 



h 




sin^^ 



— I sin 29 sin i 
\ sin 26 cos d 



\ 



1 

2 



COS 



COS 



COS 2-0 
— sin 2tlj 



sin 2ip 

cos2?/' 





sin 9 cos 
sin 6 sin 
cos 20 cos 6 sin 20 
2 V- — J) sin 20 cos ^ cos 20/ 

where {6, 0) gives the arrival direction. /i+ and h^ describe the two polarization of the 
GW in the source frame. The actual polarization in the detector depend on the position 
of the source {6, 0) and on the orientation of the source which is rotated by an angle ijj 
along the line of the sight with respect to the detector frame. We call ip the polarization 
angle since for a source at a given location a change in act only on the polarization. 

The presence of a GW in a noise-free detector leads to an output signal given by 
the vector of SQUID input current 
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where we have assumed that all Xj have the same value x ^'^^ we have defined G/, the 
sub- matrix 
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2.4- Detector sensitivity 

The detector sensitivity is given by the comparison between the output due to the noise 
and the one due to a GW. The output of the detector is given by output currents 
(one by transducer). One can combine theses outputs in different way. The ability of 
extracting the GW signal depends on the way the outputs are combined and so do the 
sensitivity. 



2.4-1 ■ Single transducer analysis In this first data analysis schema, we considered a 
rather naive way of using the different transducers: each one is taken as an independent 
experiment. This has the advantage of simplifying drastically the data analysis (which 

+ By a GW tensor we mean the perturbation of the background metric in the TT gauge 

* Note that if we describe gravitation by general relativity, only spheroidal quadrupolar modes [I = 2) 

have non zero hj 
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can be based on the one of the resonant bar without further modifications). However 
it is clear that in this way we are loosing most of the advantages of the spherical geom- 
etry. Nevertheless it is instructing to perform this analysis having in mind that it will 
probably be useful in the first phase of any sphere experiment. 

If we are looking only at the k-th transducer, its noise spectral density is given by 
the diagonal components of fl30l) 

j N 
C cjfcfe \ ^ |2c! I \ ^ I /~i |2c! 

^jfe — a — \ K:xj+2N+k,e\ >^F,e + 2^ \ ^J+2N+k,J+e\ 

£=l 1=1 
N N 

+ S^35) 

t=i e=i 
and this has to be compared with the spectral density of the k-th SQUID input 
current due to the GW. The latter is obtained as the square of the equ. fl33l) 



qGW _ 1, A 2 n2 2 
^k - ^s-^sX 




(36) 



If we know the GW arrival direction, we can compute each Ty and simulate the 
detector response. The inverse problem, computing the GW propagation direction, 
needs high signal to noise ratio (SNR) [17] . If we are interested in the detector sensitivity 
we are working at SNR=1, and therefore we have no information on the arrival direction. 
Furthermore, if we choose an arbitrary direction we can be in a case where some of the 
modes are poorly or not at all coupled to this peculiar GW. Consequently we can 
underestimate the sensitivity 0. As the spectral density is a statistical feature of the 
signal, it is then natural to perform an average on the possible arrival directions and 
over the source orientation (polarization angle). 

Averaging over the arrival direction {9, 0) and the source orientation 

(polarization angle ip) and assuming an intrinsic polarization = h, h^ = 0, we 
find [19] 

= ry^^lRlx'h' E G]^,^^,,. (37) 
1=1 

The strain sensitivity is given by the comparison = and is 
where we have defined the transfer function Tp. 



tt As an illustration: the case of a cylindrical detector. If we choose the GW arrival direction parallel 
to the detector axis, the GW is not seen and no information is given about the detector sensitivity. 
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2.4-2. Coherent analysis Rather than looking at each transducer separately we can 
perform a coherent analysis taking into account all the N signals. In this case we 
have to perform multidimensional matched filtering [17]. In the special case of TIGA 
configuration, one can use the mode channels in order to obtain directly the motion 
of the quadrupolar modes of the sphere [6] . However, on one hand, this property was 
proved for the sphere with mechanical resonators and may be affected by the tuned 
electric oscillators and/or by the perturbations induced by the suspension of the sphere 
in the Earth gravity. On the other hand, having the motion of the sphere modes is not 
enough: one has to make the deconvolution between the excitations and the damped 
motion of the modes. 

If, for a given transducer number and configuration, we know the detector response to 
an excitation, one can directly build a multidimensional matched filter [18]. One can 
therefore achieve optimal filtering (once the waveform of the GW is know or assumed to 
be know). Note that the construction of the filter can be performed using our numerical 
model, see [H] for details. 

The multidimensional matched filtering can be defined and used for any configuration 
but in general it has to be computed numerically. However the TIGA configuration 
offers the great advantage of accepting analytic solutions to the deconvolution problem. 
In the following, thanks to a result of Stevenson [17] , we do not need to know how the 
matched filtering will be performed and whether the filter is obtained in a numerical or 
analytical way. We only need to know that an optimal filter is applied to the data. In 
this case the SNR is given by [T7] 
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SNR= / a{iu)dLu/2TT (39) 
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where S is the noise matrix defined by equation (127|) . 
Using the expressions (l33l) . we get 

aiu) = lu'm'XxHK ^x)T^GjS^Tv^(^~+^ (41) 

H 

In order to evaluate this expression we face again the problem that Ty is a function 
of the arrival direction and of the source orientation. As in the single transducer analysis 
it is convenient to average the expression over the direction and the polarization. This 
gives 

2- sm{e)ded(f) d%l) (h+\ 
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assuming again a polarized GW = h, = we get for (a) 



(a) = ^u'mlR^x'h'Trin) (43) 
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The strain sensitivity (averaged on the sky and the source orientation) is then given 



3. Second part: application and simulations 

3.1. The filled sphere 

We start with the case of a filled sphere, this case is easier to compare with existing 
computations [S] and can be used as a benchmark for the hollow sphere simulations. 
We also use this simple case in order to highlight the how the single transducer 
sensitivity is modified by the presence of other transducers. With a numerical code, 
implementing f|T2|) . fll5tl2^ . fl35l) . fl38l) . and dB]), we can compute the strain sensitivity 
for a spherical resonant mass with capacitive transducers coupled to a dc-SQUlD with 
input transformer and resonances frequency of the transducer tuned in order to monitor 
the five quadrupolar modes. The parameters used in our simulation can be found in the 



3.1.1. Single transducer results In the case of a single transducer analysis, we obtain 
sensitivity curves - one for each transducer. As an example we considerate a l[m] 

radius bulk CuAl sphere with 6 transducers in TIGA configuration. At the quantum 
limit, the typical strain sensitivity of one of the transducers is plotted in Figure HI Note 
the presence of horns on the both side of the resonance. We can understand them as the 
contribution from the noises of the 5 others transducers. This explanation is confirmed 
by the sensitivity curves obtained for the same sphere but with a different number of 
transducers, see Figure O As we increase the number of transducers we also increase the 
number of noise sources, and the sensitivity of a single transducer is worse than if the 
sphere was equipped with a single transducer. Although this finding is not surprising 
and can be related to existing literature [231 121] , the details of this feature are important 
in order to calibrate the transducers. 

3.1.2. Coherent analysis results We now work out the coherent strain sensitivity for the 
same sphere as in the previous case. For 6 transducers placed into a TIGA configuration 
the sensitivity does not present horns as show in Figure [61 We also plot the sensibility 
corresponding to the best present transducer {Nphonon = 50). 

One can address the question of understanding how the sensitivity is modified if 
we change the number of transducers. As the sphere with a single transducer has 
the same pattern function as a bar, we can expect that each new transducer partially 
complete the coverage of the polarization and arrival directions until for = 5 we have 
a omnidirectional and omni-polarization detector. The sensitivity will then be improved 
by the addition of each transducer. Furthermore we will add the transducers in order 



by solving the equation (cr) = 1 for h{u!). We eventually obtain 




(44) 



Table [H 
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Sphere 


Radius 
iviass 
Quality factor 
Number of modes 
Modes frequencies 
Radial eigen-function 
Ratio effective/real radius 


Rs = l[m] 

rUg — ooKjy L\K,g\ 

Qs = 5- 10^ 
J = 5 

^ = 1049, 1053, 1057, 1061, 1065 [Hz] 
a = -2.89 
X = 0.328 


Transducer 


Mass 
Quality factor 
Resonator frequencies 
Field in the capacity 
Capacities 
Transformer primary inductance 
Transformer secondary inductance 
±iansioimei muLuai inciucLance 
Transformer resistance 


mt = 165[kg] 

Qt = 5- W 
^ = 1046, 1051, 1055, 1059, 1063, 1068 [Hz] 
i!/ — 4.0 ■ iU yV 1 Tn\ 
Ct,k = 36.9, 36.7, 36.4, 36.1, 35.8, 35.5 [nF] 
L, = l.-i[H] 
= 8 ■ 10"6[i/] 

IVl — Z.U • lU 1 -'^ 1 

Tk ~ 2.9 ■ 10-^[fi] 


SQUID 


Input inductance 
Washer inductance 

Shunt resistance 
Mutual inductance 


Li = 1.7- 10-6[iJ] 
LsQ = 80 ■ 10-''[H] 

Rsh = m 

MsQ = lO[nH] 



Table 1. Parameter values for the filled sphere at the quantum limit. If the mode 
(j) or transducer (fc) indices is absent, its means that the same value is used for all 
the indices range. The frequency splitting in LJs,j is extrapolated from the miniGRAIL 
one. The capacities Ct^k are chosen in order to tune the electric frequencies on the 
resonator ones. The transformer resistance is computed in order to get a electrical 
quality factor Qe = 5 • 10^. The temperature and the SQUID effective temperature are 
fixed to 2{)[mK]. 



to complete a TIGA configuration and therefore, due to the particular symmetry of this 
configuration, we expect the complete TIGA configuration to be even more sensitive 
than the incomplete 5-transducers one. The simulations confirm this analysis as shown 
in Figure [71 

3.2. The hollow sphere 

All the above treatment has a nice, simple generalization to the case of a hollow sphere 
[20] . While preserving all the feature of a bulk sphere such as omnidirectionality, and the 
capability to determine the source direction and wave polarization, an hollow sphere has 
several interesting peculiar properties. Its quadrupole frequencies are lower than those 
of an equally massive solid sphere, thus making the low- frequency range accessible to 
this antenna with good sensitivity. Further, as shown in [20], for an appropriate ratio 
between the inner and outer diameter, the cross section for the second quadrupole mode 
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Figure 4. Strain sensitivity at the quantum limit for one of the 6 transducers placed 
into a TIGA configuration on a l[m] radius bulk CuAl sphere. The curves are the 
relative contributions of the different noises: the thick line is the total sensitivity, the 
dashed curve is the mechanical thermal noise contribution (|20l2ip . the dashed-dotted 
line is the thermal electric noise contribution (I22I1 . the dotted line the back-action 
contribution ((23|) and the continuous curve the white noise ([24|) . 

equals that of the first, and one has the possibihty of working with a detector with the 
same (high) sensitivity at two frequencies. The main differences with the bulk sphere 
are the numerical values of the radial functions a^, that we first met in equation ([9]), 
and the one of the coupling to GW, namely x (see equation (13T|) ). We have also to 
recompute the eigen- frequencies (this can be done using the formulas from ref. [20]). 
We have applied our model to a l[m] external radius and a = 0.4 [m] internal radius 
hollow sphere, with 6 transducers in the TIGA configuration coupled to the first 
quadrupolar modes. The results of our computation are displayed in the Table [2] The 
corresponding coherent strain sensitivity at the quantum limit is show In Figure [8] Note 
that the frequency window is shifted by about 200 [Hz] with respect to the filled sphere 
and that the sensitivity loss is only by about a factor 1.5. 

3.2.1. Effect of the thickness of the hollow sphere We address now the question of the 
influence of the thickness of the sphere on its sensitivity. For a fixed external radius a 
thinner sphere is less massive, and therefore less sensitive. There is also a shift in the 
resonances to lower frequencies when the sphere is thinner, see Figure [Q] 
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Figure 5. Strain sensitivity for tlic single transducer analysis. The purpose of this 
figure is to show the effect of the noise due to the presence of other transducers. The 
dotted curve is the sensitivity for a single transducer on a l[m] radius bulk CuAl 
sphere. Then, always performing the analysis of the output of one transducer we 
add others transducers acting as noises sources: the dashed-dotted curves stand for a 
sphere with 2 transducers on the 2 first TIGA locations, the dashed curves for a set of 
4 transducers on the 4 first TIGA locations, and the continuous curve is the sensitivity 
for a transducer out of a 6 transducer TIGA configuration. Note how the presences of 
further noises sources (= the other transducers) deteriorate the sensitivity by adding 
horns. 

3.2.2. Enlarging the bandwidth by monitoring the second quadrupolar multiplet We now 
compare different designs of multi-modal resonant detectors based on a hollow sphere. 
We will always work on the base of our 0.4/1 [m] hollow sphere described above and we 
compare different transducers designs. 

Having a description of the hollow sphere with 6 transducers monitoring the first 
quadrupolar multiplet we can now simulate a multi-modal detector. As discussed in 
Section 12.1.41 we have to take into account the presence of other sphere modes. In 
Table [3l we list the properties of the modes we include in the simulation of our multi- 
modal models. Each modes multiplet (n, t) contain 2£ + 1 modes and therefore we take 
into account a total of J = 30 spheroidal modes in our simulations. Only 10 of them 
are coupled to GW. 

For a hollow sphere, we have the choice of placing the transducers on the outer or 
on the inner surface (this choice may be limited by the experimental difficulty of placing 
transducers inside the sphere). Looking only at the quadrupolar modes, and because 
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Figure 6. Strain sensitivity at the quantum limit (left) and for Nphonon = 50 (right) 
for a coherent analysis of the outputs of 6 transducers placed into a TIGA configuration 
on a l[m] radius bulk CuAl sphere. The curves are the relative contributions of the 
different noises: the thick line is the total sensitivity, the dashed curve is the mechanical 
thermal noise contribution ()20|2ip , the dashed-dotted line is the thermal electric noise 
contribution ([22|) . the dotted line the back-action contribution ((23|) and the continuous 
curve the white noise (l24l). 
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Figure 7. Coherent strain sensitivity at the quantum limit of a transducer set placed 
on a l[m] radius bulk CuAl sphere. The continuous curve is the sensitivity for a 6 
transducers TIGA configuration, the dashed curves the one for a set of 4 transducers on 
the 4 first TIGA location the dashed-dotted curves the one for a set of 2 transducers on 
the 2 first TIGA location and the dotted curve is the sensitivity for a single transducer. 
(Note that the dotted curve is the same as in Figure [5]) 
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Sphere 


iXcxU-lUo 

Internal radius 
Mass 
Modes frequencies 
Radial eigen-function 

J. VjGi V±\J x^jLLkj^^ u 1 V / L v3 ( X i L ( X \ LI. U.O 


JTLs — 

a = 0.4[m] 
rus = 30974[A;^] 
^ = 861, 865, 869, 873, 877 [Hz] 
a = -2.73 
V — n 322 


Transducer 


Mass 

Resonator frequencies 
Field in the capacity 

Capacities 
Transformer resistance 


nit = I55[kg] 
^ = 859, 863, 867, 871, 875, 879 [Hz] 
E = 4.5- W[V/m] 
Ct,k = 55, 54, 54, 53, 53, 52 [nF] 
Tk ~ 2.4- 10-^[fi] 



Table 2. Parameters value for the hollow sphere with a = 0.4[m] internal radius at 
the quantum limit. The transducer being build on the same model as for the filled 
sphere we put in this table only the quantities which differ from the filled case. The 
temperature and the SQUID effective temperature arc fixed to 2Q[mK]. The other 
values can be found into the Table [TJ 




700 750 800 850 900 950 1000 1050 



frequency [Hz] 

Figure 8. Coherent strain sensitivity at the quantum limit for a set of 6 transducers 
placed in a TIGA configuration on a 0.4[m]/l[m] internal/external radius CuAl sphere. 
The curves are the relative contributions of the different noises (see FigurelH). The gray 
curve is the sensibility corresponding to the best present transducer {Nphonon ~ 50). 
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Figure 9. Coherent strain sensitivities at the quantum limit for a TIG A configuration 
on a 1 [m] external radius hollow CuAl sphere for different values of the internal radius 
a. The dotted curve is for a ~ 0.8[m] (ie. a shell of thickness 0.2[to]) the successive 
curves are for a decreasing by step of 0.2[m] until the black thick line which stand for 
the filled sphere. 





fn,l 


ttn.K-Rs) 


an,i(o) 


Xn,l 


(1,2) 


869 [Hz] 


-2.73 


-2.85 


0.319 


(1,3) 


1421 [Hz] 


-0.488 


17.1 





(1,1) 


1558 [Hz] 


-1.41 


-4.02 





(1,0) 


1758 [Hz] 


-1.29 


-1.06 





(1,4) 


1937[Hz] 


1.51 


21.1 





(2,2) 


1970 [Hz] 


0.276 


-3.27 


-0.148 


other 


>2.4[kHz] 









Table 3. Frequency, radial eigen-function (at the outer and inner surface) and coupling 
to GW for the first modes of a hollow sphere with external radius Rs = l[m] and 
internal radius a = 0.4[m]. All these quantities are computed from formulas found 
in [2D]. 
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the radial displacement of the second mode is maximum at the inner side of the sphere 
[20] . one could think of improving the sensitivity by mounting the classical capacitive 
transducers considered here inside the sphere. However, as show in Table [3l the higher 
multipolar modes have very large a's on the inner surface. These modes will couple 
to the transducers main resonances and their thermal noise will eventually limit the 
detector sensitivity. This effect is analysed in detail below. 

A first approach to design a multi-modal spherical detector is to add a second set 
of 6 transducers in TIG A configuration coupled to the second quadrupolar multiplet. 
The characteristic of the second set of transducer are listed in Table HI The parameters 
are optimized to get the best strain sensitivity and the largest bandwidth following the 
detailed analysis for a bulk sphere reported in fi2[ [13] . We notice that the optimal 
mass of the resonator coupled to the second quadrupolar modes is smaller that the first 
resonator. This is due to the fact that the effective mass of the higher quadrupolar 
modes is smaller and the higher frequencies of the modes requires a smaller transducer 
mass to have the same electro-mechanical coupling. In our simulation we can choose to 
place the transducers inside or outside the sphere. On Figure [TOl we show the difference 
in the strain sensitivity for a hollow sphere with 12 transducers when the thermal noise 
contribution of all the modes between the first and second quadrupolar multiplet is 
included. Both the configurations with the second set of transducers placed outside 
and inside the sphere are considered. The thermal noise of the modes not coupled to 
GW reduces slightly the bandwidth in the former case, but has dramatic effect on the 
sensitivity in the latter case. 



Mass 


m,k=7..12 ~ |"^t,fc=1..6 = 31[%] 


Resonator frequencies 


^ = 1947, 1956, 1965, 1974, 1983, 1992 [Hz] 


Field in the capacity 


E = A.5- 10''[V/m] 


Capacities 


Ct,k ~ 10 [nF] 


Transformer resistance 


Tk ~ 5.4 ■ 10-^[fi] 



Table 4. Parameter values for the second set of transducers. For comments, see 
Table [T] and for the sphere parameters see Table [2] and [3] . 

From Figure [10] one can conclude that placing the transducers inside the sphere 
will not lead to an improvement of the sensitivity. We observe that a 12 transducers 
in 2 TIGA configuration on the outer surface of a 0.4/1 [m] hollow sphere can lead to 
a sensitivity curve with 2 windows below 3 ■ 10~^^[-frz~^/^] with bandwidth of 175 and 
52 [Hz]. 

Having in mind that the presence of 12 transducers can represent an experimental 
problem (even on the outer surface of the sphere), we now consider the possibility of 
reducing the number of transducers by using the double-mode transducer described in 
Section 12.1.31 We then perform simulations of a multi-modal spherical GW detector 
with only 6 double-mode transducers in a TIGA configuration. 
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Figure 10. Coherent strain sensitivity at the quantum limit for two sets of 6 
transducers placed in TIGA configuration on a 0.4[m]/l[m] internal/external radius 
CuAl sphere. Each transducer set is coupled to the first and second quadrupolar 
modes, respectively. The curves show the relative contributions of the different noises: 
the black thick line is the total sensitivity, the dashed curve is the total mechanical 
thermal noise contribution (|20|2ip , the dashed-dotted line is the thermal electric noise 
contribution ([22|l . the dotted line the back-action contribution ([23|l and the continuous 
curve the white noise ((24|) . The left panels are the result of the simulations with only 
the quadrupolar sphere modes, in the right panel we have add the other modes with 
frequency smaller than 2.5[kHz] The upper plots are for transducers on the outer 
surface and the lower ones are for transducer on the inner surface. On the upper right 
plot the gray curve is the sensitivity for Nphonon — 50 



In figure [TT] we show the sensitivity of a hollow sphere equipped with 6 double-mode 
transducers. The transducers parameters, reported in Table [5], are chosen to get the best 
sensitivity and the largest bandwidth at the two principal quadrupolar modes multiplets 
around 800 and 1900 Hz. The choice of the parameters in Table [5] is the result of an 
iterative process based on the detailed optimization procedure described in [T2l [T3] for 
a bulk sphere. A more detailed technical analysis is necessary to fully optimized the 
detector described here. Such an analysis is beyond the scope of this paper. 

The sensitivity around the two most sensitive spheroidal modes is comparable with 
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' l,k 'J 




r2,fc ~ 5.4 ■ 10-4[fi] 


Input inductance 


Li = 1.6- 10-6[i/] 


Washer inductance 




Shunt resistance 


Rsh = m] 


Mutual inductance 


MsQ = lQ[nH] 



Table 5. Parameters value for the set of double modes transducers. For comments 
see Table [1] and for the sphere parameters see Table [2] and [3l 



the one obtained with a 12 single-mode transducers configuration. Thanks to the use 
of a single SQUID to amplify the signal from both the spheroidal modes families, the 
resulting back-action noise contribution of the SQUID amplifiers is reduce in the double- 
transducer read-out scheme. This leads to a larger bandwidth at the second spheroidal 
mode frequencies. Between the two modes families the sensitivity is reduced due to a 
low mechanical coupling between the resonators and the sphere modes, and is limited 
by the additive SQUID white noise. A strain sensitive level of about 10~^^ Hz~^/^ 
and a total bandwidth of about 600 Hz could be achieved in this way around 800 and 
2000 Hz. At the sensitivity of the currently best performing resonant bar antenna [1] , a 
multimodal hollow spherical detector could reach a bandwidth of about IkHz. 

3.3. Concluding remarks 

We developed a mathematical framework, which can describe the indirect coupling of the 
transducers through the modes of the sphere. We have shown how to build models giving 
the strain sensitivity of a GW spherical detector and presented different configurations 
of a multi-mode spherical resonant detector. The effect of the thermal noise from 
higher frequency modes of the sphere, not coupled to GW, are also included in the 
model. We fully analyzed the strain sensitivity of an hollow CuAl sphere of 0.4[m]/l[m] 
internal/external radius. Such an antenna, equipped with 12 'standard' transducers or 
6 double-mode transducers, when optimally tuned, may display a sensitivity curve, at 
the quantum limit, below 3 ■ 10"^^ [/J 2;^^/^] with a total bandwidth up to 600 Hz around 



On the Sensitivity of a Hollow Sphere as a Multi-modal Resonant Gravitational Wave Detector^?) 



10 



-18 



N 10' 




10 



-26 



500 



1000 



1500 

frequency [Hz] 



2000 



Figure 11. Coherent strain sensitivity at the quantum limit for a set of 6 double- 
mode transducers placed in TIGA configuration on a 0.4[m]/l[m] internal/external 
radius CuAl sphere. The curves show the relative contributions of the different noises: 
the black thick line is the total sensitivity, the dashed curve is the total mechanical 
thermal noise contribution (|20|2ip , the dashed-dotted line is the thermal electric noise 
contribution ([22|) . the dotted line the back-action contribution ((23|) and the continuous 
curve the white noise (|24p . All the spheroidal modes with frequency smaller than 
2. 5 [kHz] are included in the simulation. 



the first two quadrupolar multiplets at about 800 and 2000 Hz. At the sensitivity of 
the currently best performing resonant bar antenna [1], a multimodal hollow spherical 
detector could reach a bandwidth of about IkHz. Such a detector could have peak 
sensitivity comparable with the first generation of interferometers, but would be also 
able to determine the GW arrival direction. With an appropriate resizing of the sphere, 
the second quadrupolar modes can be shifted to the frequency region where existing 
small spherical detectors are sensitive. This would open the possibility of coincidence 
search between several spherical detectors and the DUAL detector [21], building in this 
way the base for a powerful omnidirectional gravitational wave observatory. 
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